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Abstract 



• With the generalized factorization approximation, we calculate the branch- 

> 

ing ratios and CP asymmetries in B meson decays into two charmless 
pseudoscalar mesons. We give a new estimation of the matrix elements of 

o ' 

CT\ ■ (S + P)(S — P) current products with the PQCD method instead of equation 

ON 

r~| ■ of motion. We find that our results are comparatively smaller than those in 

a: 



the literature. 
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I. INTRODUCTION 



Penguin diagrams play an important role in charmless B decays and direct CP violation. 
They can provide not only the necessary different loop effects of internal u and c quarks 
|U , but also dominate the branching ratios of many modes of charmless B decays, such as 
B tvK, B -» KK and B -> Krj' . 

As we know, the standard theoretical framework of studying non-leptonic B decays is 
based on the effective Hamilton approach and the factorization approximation. The effective 
Hamiltonian is expressed by a sum of the products of a series of Wilson coefficients and four- 
quark operators. Unfortunately, we have many difficulties in calculating matrix elements 
of the four-quark operators directly in exclusive non-leptonic B decays, such as B to two 
pseudoscalar meson processes. So we have to use the factorization assumption, usually the 
BSW model ||. Then the mesonic matrix elements are factorized into the product of two 
matrix elements of single currents, governed by decay constants and form factors. However 
in the BSW model, the factorization involves the contributions of Fiertz transformations 
of the four-quark operators. Using the Fiertz rearrangement, one can find that the current- 
current product (S + P)(S — P) matrix elements should be taken into account. The general 
method to deal with (S+P)(S—P) matrix elements is to transform them into (V— A)(V— A) 
matrix elements by using equation of motion. Then one can find that (S + P)(S — P) matrix 
elements are very sensitive to the masses of light quarks. But the current masses of light 
quarks are not determined precisely. Obviously, it brings large uncertainty for estimating 
the (S + P)(S — P) matrix elements and the branching ratios of charmless B decays. As 
pointed by A. Ali, G. Kramer and CD. Lii H, varying the light quarks masses by ±20% 
yields variation of up to ±25% in some selected decay modes(such as — > K ± rj and 
B° — > K°7] .) In this work, instead of using equation of motion we will try to apply the 
PQCD method to recalculate the ratio of the (S+P)(S — P) to (V-A)(V-A) matrix elements 
at leading twist approximation, which is not sensitive to the masses of light quarks. We think 
that it might have less uncertainties than the results obtained by using the quark equation of 
motion. So, it is necessary to recalculate the branching ratios for B meson charmless decays 
by using PQCD method and compare with those by using equation of motion. On the other 
hand, nonfactorizable effects in charmless B decays can not be neglected. To compensate it 
for this, the general approach is to replace the number of corlors N c by a phenomenological 
color parameter Nj:**. We will discuss the differences while iV'" equal 2, 3, oo respectly. 

This work is organized as follows: Section 2 gives the framework of calculation including 
the effective Hamiltonian and fatorization approximation. Section 3 is devoted to the PQCD 
method to estimate the ratio of the (S ± P)(S — P) to (V — A)(V — A) matrix elements. 
In Section 4, we calculate the branching ratios of charmless B decays into two pseudoscalar 
mesons and their CP asymmetries using the method mentioned above. We also give some 
discussions of the numerical results. Section 5 is for the concluding remarks. 

II. CALCULATIONAL FRAMEWORK 

The \AB\ = 1 effective Hamiltionian is 
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£ v q CMQ'M + CMQ q M + £ cmqm 



q=u,c 



k=3 



+ h.c, 



where w g = V^V^for 6 — > transition) or t> g = V^V^* (for 6 — > s transition) and Cj(/x) are 
Wilson coefficients which have been evaluated to next-to-leading order approximation. In 
the Eq.(l), the four-quark operators Qi are given by 



Qi 

Qi 

Q9 



(u a bf3)v-A(qf3U a )v-A 

{u a b a ) v _ A (qpU(]) v _ A 

(?A)v-AE(?Jgi)v-i 



(q a b a ) V - A Y,(q'pqp)v+A 

q' 



QI 

Qe 



{c a b/3)v-A{qj3Ca)v-A 
(c Q b a ) V -A(q(3C/3)v~A 

{qpb a ) V - A Y.{q a qp)v-A 



(qf3b a ) V -AE(q a qf3)v+A 

q' 

liqMv-AYlZq'iqpqplv+A Qs = l{qpb a ) V - A Y^e^(q a q'p) v +A 



|(3o6h)v-aEv(9J?J)v-a Qw = \(qpb a )v-AY<eq>(q a qi3, 



V-A 



(2) 



with Q\ and Q\ being the tree operators, Q 3 — Q G the QCD penguin operators and Qi — Q w 
the electroweak penguin operators. With the renormalization group method, we can evolove 
the renormalization scheme independent Wilson coefficients C*j(/i) from the scale \x = mw 
to n — hxlGeV ~ ms, which are 

C_x = -0.313 C_2 = 1.150 C 3 = 0.017 C_ A = -0.037 

C 5 = 0.010 C 6 = -0.046 C 7 = -0.001a em C 8 = 0.049a em (3) 

C 9 = -1.321a em Cio = 0.267a em . 

So we can express the physical amplitude as follows 

<Q» ■ C(a*)) = (Q T )o • C(a*) (4) 

where (Q T )o denote the tree level matrix elements and 

1 — Oi 2 — O2 O3 — O3 jj- 4 — O4 + -Tg 

C 5 = C5 — ^ Cg = + P s Cy = C 7 + P e Cg = C 8 (5) 
C 9 = Cg + P e C{ = Cio, 
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P *=^M (f GK,?V)) 



^(3^) + ^)) (y-GK?V)) (6) 



G(m q , q , (i) — —4 dx x(l — x)ln 



' 2 ^/i ^„2' 







— X{1 — X)q 



As noted in the introduction, we have to calculate the matrix elements of the four- 
quark operators by using the factorization assumption. Here we apply the BSW model 0. 
However, nonfactorizable effects are not negligible in the process of B to two light mesons. 
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We will use the simpliest approach to compensate it by using only one color parameter N*ff, 
even if there is no reason why using only one single parameter N^f to explain the branching 
ratios of all kind of different modes. 

For illustration, we give the amplitude of B~ — > K~r]' as an example: 

(K- v '\n eff \B-) 

G 



-7= S v q [(a 1 5 gu + a 3 + a 9 )Mf u / 

V ^ q=u,c 



+(a 2 5 qu + 2a 4 - 2a 6 - | + >)M^" 
+ (a 4 -a 6 + a 3 + y + y- -yjACL 



+ (a 2 5 qc ~a 8 + a w )M^ 

+ (-2a 5 - 2a 7 )XfJ + (-2a 6 + a 7 )X^~} , (?) 



where 



«2i-l — C 2i + J^ff > a 2* ~~ ^2i-l + , re 2 ff 5 (8) 

N c iVc 



and 



<««, = (^l(^2) y _ A |0)(Pl( ? - 3 6) v _ A |S) 



-PP' 



<L 3 = (^l(^2) 5+ p|0)(P'|( g - 3 6) 5 _ P | J B). (9) 



We will use the following parameterization for decay constants and form factors: 
{0\V li -A li \P(q))=if P q li 

(P 2 (q 2 )\V, - A^fa)) = F Pl ^ P2 (g!)g + , + P_ p ^ P2 (g!)g-„ (10) 

where g± = q± ±q 2 , and we use the monopole dominance assumption for the q 2 _ dependence 
of the form factors: 

F p ^(q 2 )^ - mi ~ m2 F^(q 2 ). (11) 
yH 1 m l +m 2 + yH 1 v 1 

Then we can obtain 



^PP' . = "^P JufPP' 

919293 (( mi + m2 )( m3 _ mfe )) 9!9293' 
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where Eq.(13) is derived from the equation of motion and presents the mass of the light 
quark t& respectly(z = 1, 2, 3). 

Calculation in this framework have been discussed in detail in some papers HQ involving 
the branching ratios and CP asymmetries in Non-leptonic charmless 2-body B decays. In 
these papers the uncertainties resulting from the renormalization scale dependence, non- 
factorizable contributions and the input parameters (a s , quark masses and form factors) have 
been worked out. Further penguin effects and the strong sensitivity of the CP asymmetries to 
the CKM parameters (p, 77) have been discussed there. In these uncertainties, the uncertainty 
of light quark masses is mainly showed in the part of (S + P)(S — P) matrix elements. 
Sometimes they are dominant terms in some modes of charmless B decays, such as B~ — > 
K~i]', in which the term X%J£ is enhanced by the factor Then it motivates us to give 
a new estimation of the matrix elements X^ 2qz to cancel the uncertainty of the light quark 
masses. 



Brodsky et al. |J has pointed out that the factorization formula of PQCD can be applied 
to the exclusive B decays into light mesons for the large momentum transfers. One can write 
the amplitude as a convolution of a hard-scattering quark-gluon amplitudes (j)(x, Q 2 ) which 
describe the fractional longitudinal momentum distribution of the quark and antiquark in 
each meson. An important feature of this formalism is that, at high momentum transfer, 
long-range final state interactions between the outgoing hardrons can be neglected. In the 
case of non-leptonic weak decays the mass squared of the heavy meson rn H establishes the 
relevant momentum scale Q 2 ~ rn 2 H , so that for a sufficiently massive initial state the decay 
amplitude is of the order of a s (Q 2 ), even without including loop corrections to the weak 
hamiltonian. The dominant contribution is controlled by single gluon exchange. 

We intend to apply PQCD method to estimate those hadronic matrix elements such as 
(V — A)(V — A) and (S + P)(S — P) at the leading twist approximation. The wave function 
of B meson and flavor SU(3) singlet or octet pseudoscalar mesons are taken as: 



where Iq is an identity in color space. In QCD, the integration of the distribution amplitude 
is related to the meson decay constant 



III. PQCD METHOD 



= J-ICL 

= 1 kt 

V2 V3 



4p{y)(i + mph5 , 



(14) 




(15) 



Then we can write down the amplitude of Fig.l as 




{P\qi7bQ»\0) pqcd = 3 x 





dy4>p(y)Tr [75 (i + m P )7 5 ] = f P m P 



(16) 
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In a consistent way, we can use perturbative QCD to estimate the matrix elements like 
(PlqijfjblB) and (P\q~ib\B) (Fig.2, Fig. 3), where q\ denotes light quark field operator and we 
have neglected the fermi motion of quarks, while the gluons in the Fig. 2, 3 are hard because 

k 2 = ( X p - (1 - y)q) 2 ~ -x(l - y)m 2 B ~ lGeV 2 (17) 

(here we using mean values (y) ~ |, (x) ~ e B , with e B ~ 0.05 — 0.1 and x 2 << 1) so, we can 
neglect the 0(x 2 m B ) term and use perturbative QCD method to calculate the amplitude. 
It turns out to be 

2 f 

{P\qa^b\B) pqcd = -^9 2 J dxdy<j> B (x)<l>p(y) (18) 
jTr [75 + m P )Ytfan(ft + m B)lslu] Tr [75 + m P )"j ^ + m h )Y(i> + ™b)1s1v} 1 



I k2p i fc 2 (P 6 2 -mg) j ' 

In order to get quantitative estimation, we take the wave functions as P,[7| 

MX) = 276 5{X ~ ' MV) = ] H fMl ~ V) ' (19) 

(here e B is the peaking position of the B-meson wave function. Typically (e B ) ~ mB m™ b ) 
We get 

(P|ft7 M 6|S) pqcd = F^ p {Q 2 ){p + g)„ + F^ p (Q 2 )(p - g)„ (20) 



where 



^B-,p (n .2\ 8na s , , I m P m B 

+ (g) -~~ /p/B \ 4< 

m b (m P - 2m B ) + y(m| - 2m P m B ) \ 
£Bm B (ym B - m b ) J 



2 4 

2m b - 4m P - y(m B - 2m P ) 1 

*» wvmi - — I' (22) 

Here, Q 2 = (p — g) 2 .So we can obtain the matrix element . We can also get the matrix 

element (P\q~ib\B) as 



2 C 

(P\sb\B) pQ CD = --g 2 J dxdycp B (x)<pp(y) . ^ 2p2 

Tr [jsjj + m P )(P b + m b )Y{ff + m B )^ v ] \ 
k 2 (P 2 -m 2 ) } 

8tt j -2m P (l - 2e B ) + e B m B )m B + nip 
- - Y a s f B f P j ^3 

m B (m P + m b ) — 2m 2 R — 4m b m P + ympm B ] 

dyv w^i-mi) r (23) 
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In the literature 0, as an example the authors calculated the numerical results of the 
matrix element (K~\s~f IJ ,b\B~) by above framework, where they applied a s ~ 0.38, f B = 
200MeV and fx = 160MeV. One can find their results are sensitive to the values of 
parameters e B and m b , and seem small compared with the BSW result. We also compute 
the matrix element (ir~\s , -f fl b\B^) and list the numerical results in Table 1. Where we take 
a s = 0.38, f B = 0.2GeV and / w = O.lZGeV. 

We can see that the results are very sensitive to the values of parameter e B and nib, and 



smaller than the BSW result which is about 0.29 [14|. As mentioned in the Ref. |J, the 
PQCD results are comparatively small in many cases. 
But the ratio 

X pp " 

<t? _ <?i<?2g3 PQCD ( 0A \ 
M pp ' 1 ' 

m qiqiqz PQCD 

is insensitive to the parameters e B and m\>. So it is more reliable because of cancelation of 
the main uncertainties. We list our computation in the Table 2. 
The ratio by using equation of motion is 

X K°n- 2_ 

n= ±^—= Hhf ---0.30, (25) 

M ™ 3 (m s + m d )(m d -m b ) 

and it is about one order of magnitude larger than the PQCD estimation. As mentioned 
in our introduction, the matrix elements of (S + P)(S — P) four quarks operator are very 
important in some decay modes of B mesons, like B to rj' and other mesons. So it is necessary 
to recalculate the branching ratios and CP asymmetries for 2-body charmless B decays by 
using the PQCD method instead of the equation of motion. 



IV. BRANCHING RATIOS AND CP ASYMMETRIES 

In the B rest frame, the two body decay width is 

1 i/oo/iet imia \P\ 



r(£? _> P pi) = —\(pp'\ Heff \B)\ 2 ^, (26) 

where 



8vr n m 2 B 



= [{rn? B - (m P + m P ,) 2 )(m 2 B - {m P - m P ,) 2 )} 2 

is the magnitude of the momentum of the particle P or P'. The corresponding branching 
ratio is given by 

Y(B — >■ PP') 

B BR {B - PP') = (28) 

t tot 

The direct CP asymmetry Acp for B meson decays into PP' is defined as 

= YiB ^ PP') - T(B -> PP) 
CP T(B — > PP') + V(B — >■ PP') ' y } 
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In our numerical calculation, we use the Wolfstein parameterization for the CKM matrix 



V ( 



CKM 



1 - 



2 



-A 



A 
1 - 



2 



\ 3 A(l-p-ir]) —\ 2 A 



X 3 A(p — irj) 
X 2 A 
1 



+ 0(A 4 ) 



and we take 



A = 0.823 ± 0.033, A = 0.220, 



P 



0.160, 



n 



0.336. 



(30) 



(31) 



Otherwise we take all parameters such as meson decay constants and form factors needed 
in our calculation as follows g[T3J: / w = O.UGeV, f K = 0.160GW, f u , = f d , = 

-0.0063GeV, 



0.049GW, f% = 
fc = -0.0024GW, and /: 
^'(0) 



0.092GeV, 



O.UGeV, /* = -0A05GeV, f°, 



(0) 



0.29, 



2x0.282 



V6 



'(0) 



0.307 



F_ 
♦"(0) 



(0) = 

0.335 



0.32, F. 



B,, 



+ " (0) 



0.254 
V6 ' 



wave functions of rj and as 



IV) 



\uu) + \dd)+2\ss) 
\uu)+\dd) — \ss) 

73 • 



Here we apply the flavor 



(32) 



We give the numerical results of the branching ratios and CP asymmetries for B charm- 
less decays in Table 3,4,5. As a comparison, the results by using the equation of motion 
are also listed in the tables where we take m u = 5MeV, = lOMeV, m s = 150MeV and 
mj = 5.0GeV. In the calculation, we have neglected the contributions of W-annihilation, 
W-exchange and space-like penguin diagrams. 

From the tables, we can see the following features: 

(i) For most of charmless B decays, the contributions of penguin diagrams are important. 

(ii) Comparing the results of the PQCD method with those by using equation of motion, 
one can find large difference between them. In the modes of B — > ttK and B — > KK, the 
branching ratios predicted by using equation of motion is larger than those of the PQCD 
method by about a factor 2. While final state involving rj , the factor would be more large. 
Obviously, CP asymmetries are also affected by these difference. In our computation, we 
find that the ratio of to M^? 2q3 predicted by PQCD method is not of m 2 P dependence 
like the estimation by use of equation of motion. So while mp is large, the distinguishes 
between two method are more obvious. 

(iii) In many decay modes, the branching ratios are sensitive to the color parameter N%ff, 
such as Bj — > 7r°7r°, 7777(77 77 ) and 7777 . Otherwise, the value of N%" affects CP asymmetries 
more largely than branching ratios in some modes, for example, B® — > K°i] , which CP 
asymmetry rangs from 60.6% to —50.1% for N^'* ranging form 2 to 00. It is because that 
dj are sensitive to N%* f which gives the different strong phases. 

(iv) Our results are smaller than those in some literatures P,p!3|. In some decay modes 
like B — > Krf, our results are one order of magnitude smaller than the results of the 
experiments |[L0|| . Because we did not consider the contributions of other mechanisms, such 
as b — > sg* — > sr] via QCD anormaly [HJ, b — > sgg — > sr] f [p| etc.. In the Ref. Rp^, the 



authors gave the numerical results involved the contributions of the new mechanisms, which 
fit the experiments very well. 
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V. CONCLUDING REMARKS 



In this paper, we recalculate the decays of B to two charmless pseudoscalar mesons 
with conventional method (the standard effective weak Hamiltonian and the BSW model). 
Instead of using equations of motion, we use an alternative method to estimate the hadronic 
matrix elements (S + P)(S — P) and obtain comparatively smaller results. In some modes, 
which are penguin dominant, such as B — > nK, the branching ratios that we predicted seem 
to be a little bit smaller than the lower limits of the experiments of CLEO |IOfl . But they 
are derived in the factorization approach, many mechanisms are not considered in this work 
such as final state interactions. Especially in the modes of B — > ttK or KK, FSI could 
yield dominant contribution to the decay width ||. So more uncertainties in non-leptonic 
charmless B decays need us to study in the future. 
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TABLES 





e B = 0.05 


e B = 0.06 


e B = 0.07 


e B = 0.08 


mb = 5.0GeF 


0.21 


0.17 


0.14 


0.12 


mb = 4.9GeV 


0.26 


0.21 


0.18 


0.15 


mb = 4.8GeV 


0.19 


0.15 


0.13 


0.11 



TABLE I. The PQCD estimations about the element (7r _ |s7 M &|.B-). 



Y K°n- 
-p add PQCD 

m sdd PQCD 


e B = 0.05 


e B = 0.06 


e B = 0.07 


e B = 0.08 


mb = 5.0GeV 


-0.025 


-0.025 


-0.026 


-0.026 


mb = 4.9GeV 


-0.025 


-0.026 


-0.026 


-0.026 


mb = 4.8GeV 


-0.026 


-0.026 


-0.027 


-0.027 



M 



TABLE II. The PQCD estimations about the ratio of the matrix element X^cTpqcd *° 

K°TT- 

sdd PQCD- 
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TABLE III. Branching ratio in 10 -5 ,and CP asymmetries in %. 'QCD' and 'EW present the 
QCD penguin and EW penguin effects respectively, and 'DIRAC presents the results with the 
equation of motion. 



N%ff = 2 






Branchin 


g Ratio 








CP Asymmetry 


Decay Mode 


TR 




QCD 


EW 




DIRAC 




QCD 


EW 


DIRAC 


B~ - 7r°7r- 


0.53 




0.53 


0.53 




0.54 




0.03 


1.9 


3.5 


B u *~V 


0.14 




0.19 


0.19 




1.54 




16.4 


0.05 


0.2 


B u *~V 


0.42 




0.57 


0.57 




1.48 




16.0 


0.1 


0.6 


B- -+ K°K- 






0.032 


0.032 




0.065 




12.9 


13.3 


12.3 


B~ -+ tt°K- 


0.038 




0.090 


0.19 




0.405 




-39.0 


-23.3 


-14.8 


B~ - K-r,' 


0.17 




0.42 


0.38 




1.38 




-11.5 


-16.7 


-7.8 


B- - 


0.024 




0.063 


0.038 




0.025 




17.8 


-7.4 


-10.4 


S" - K°7T- 






0.35 


0.33 




0.72 




-0.3 


-0.3 


-0.3 


B° d - vr+vr- 


0.64 




0.76 


0.83 




0.77 




13.2 


13.2 


18.1 


50 -► TfV 


0.012 




7.9 x 10~ 3 


6.8 x 10- 


-3 


6.94 x 10 


-3 


-42.7 


-46.8 


-51.5 


32 - A' 


1.6 x 10" 


-3 


7.0 x 10~ 3 


5.0 x 10" 


-3 


5.0 x 10" 


3 


-27.8 


-36.0 


-36.0 


Bj> - A 


3.1 x 10" 


-4 


0.019 


0.019 




0.34 




0.8 


0.8 


8.6 




1.1 x 10" 


-3 


1.2 x 10~ 3 


1.1 x 10" 


-3 


8.63 x 10" 


-3 


7.5 


5.6 


12.2 


£° W 


9.4 x 10" 


-3 


7.6 x 10- 3 


8.1 x 10" 


-3 


8.47 x 10" 


-3 


2.2 


4.8 


35.5 


5° -► 777/' 


0.056 




0.058 


0.056 




0.094 




9.3 


9.3 


95.7 


5° -► 






0.042 


0.041 




0.065 




12.8 


12.6 


12.3 


B° d - 


0.048 




0.16 


0.23 




0.46 




-31.7 


-28.5 


-18.2 


5° -► ir°K° 


1.3 x 10" 


-3 


0.22 


0.14 




0.38 




5.3 


8.7 


3.9 


B° d -► i?°r/ 


0.017 




0.42 


0.36 




1.47 




-2.9 


-10.8 


-5.4 




2.4 x 10" 


-3 


0.013 


3.6 x 10" 


-3 


2.8 x 10- 


3 


13.7 


-45.2 


2.3 




0.68 




0.81 


0.81 




0.89 




13.1 


13.1 


18.1 




0.020 




0.012 


0.012 




0.012 




-40.9 


-44.7 


-45.8 


4° - kS 


8.2 x 10" 


-3 


0.063 


0.060 




1.16 




62.2 


60.6 


19.5 


B° s K°r/ 


0.022 




0.019 


0.019 




0.34 




-8.7 


-4.9 


34.9 


50 ^ K o^o 






0.35 


0.33 




0.75 




-0.3 


-0.3 


-0.3 








0.31 


0.35 




0.75 




-0.3 


-0.3 


-0.3 


4° - A' 


6.1 x 10" 


-4 


6.1 x 10~ 4 


4.3 x 10" 


-3 


4.3 x 10- 


•3 











5° -► TT°T] 


4.6 x 10" 


-4 


4.6 x 10" 4 


3.3 x 10" 


-3 


3.3 x 10- 


3 











B»^i]r) 


7.4 x 10" 


-3 


0.11 


0.11 




0.41 




-3.7 


-3.8 


-1.7 


£° -► 777/ 


8.2 x 10" 


-4 


0.078 


0.078 




0.11 




4.8 


4.6 


3.2 


B° - w 


3.1 x 10" 


-4 


0.33 


0.33 




0.82 




1.4 


1.2 


0.4 
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TABLE IV. Branching ratio in 10 -5 ,and CP asymmetries in %. 'QCD' and 'EW present the 
QCD penguin and EW penguin effects respectively, and 'DIRAC presents the results with the 
equation of motion. 



AT e ff — 3 






Branchinj 


j rxaiio 








CP Asymmetry 


uecay ivioae 


1 It 




QCD 




th VV 








QCD 


EW 


DIRAC 


R — , _o_ — 


U.4Z 




0.42 




n a o 

U.4Z 




n a q 
U.4o 




0.04 


1.9 


3.8 


B~ -> 7r~r/ 


n 1 n 




0.15 




U.lO 




1.0b 




18.9 


0.04 


0.2 




U.ol 




0.46 




U.4D 




i.oy 




19.2 


0.2 


0.7 


R — , 7^0 XT' — 






0.037 




U.UoD 




U.U ( 




12.8 


12.9 


12.1 


R — , _0 TS — 


n nqi 
U.Uol 




0.12 




U.ZZ 




n ar 

U.40 




-32.3 


-20.4 


-12.9 


R — , JS—J 


o.y x ru 


-3 


0.42 




U.oo 




i /i/i 

1.44 




-11.0 


-15.9 


-7.2 


B u ^ K V 


n no i 
U.Uzl 




0.067 




n riQn 
U.Uoy 




n noi 
U.Uzl 




15.4 


-11.1 


19.8 


R— , TS0„ — 

B u —> K it 






0.41 




n a n 
U.4U 




U.o4 




-0.3 


-0.3 


-0.3 


r0 _^ ^+^- 

JD d — > 7T 7T 


U. f 1 




0.85 




n or; 




n no 

u.yz 




13.3 


13.3 


18.3 




0.0 X 1U 


-4 


2.8 x 10" 


-3 


1.4 X 1U 


-3 


O.O X IU 


3 


-31.8 


-50.3 


-30.3 


r0 _^ _0 ' 


1 o w i n- 
l.z X 1U 


-4 


9.3 x 10" 


-3 


b.y x 1U 


-3 


o.y x iu 


3 


3.2 


2.3 


2.3 


B d ^7T 7] 


Z.Z X 1U 


-5 


0.023 




n noQ 
U.Uzo 




U.OO 




10.4 


10.8 


10.8 


R _^ ' ' 


0.1 X 1U 


-5 


1.1 x 10" 


-4 


(.0 X 1U 


-5 


U.Ulo 




28.3 


0.3 


1.1 


B d -> W 


0.0 X IU 


-4 


4.1 x 10" 


-4 


/.y x iu 


-4 


U.U1 ( 




13.2 


61.7 


12.1 


r0 , 


4.U X 1U 


-3 


4.7 x 10" 


-3 


/I o v i n~ 
4.Z X IU 


-3 


U.U / o 




42.2 


42.8 


21.0 


B d ^ K K 






0.049 




U.U40 




U.U ( D 




12.4 


12.5 


12.1 


rO , _+ j^— 

— > 7T A 


U.UOo 




0.23 




n o^ 
U.ZO 




U.OU 




-30.4 


-29.1 


-18.7 


— > JT A 


y. / x iu 


-5 


0.24 




U.lO 




n /i o 
U.4Z 




1.1 


1 8 

J. . O 


7 




1.2 x 10" 


-3 


0.43 




0.38 




1.56 




-1.1 


-7.2 


-3.7 


B° d -► K°r/ 


1.7 x 10" 


-4 


0.015 




2.8 x 10" 


-3 


4.9 x 10- 


3 


2.8 


-67.9 


-78.5 




0.75 




0.90 




0.90 




0.98 




13.3 


13.3 


18.3 




1.4 x 10" 


-3 


4.3 x 10" 


-3 


2.0 x 10" 


-3 


2.2 x 10 _ 


3 


-32.6 


-52.3 


-49.6 


4° - X°t/ 


5.9 x 10" 


-4 


0.043 




0.041 




1.2 




30.9 


29.0 


13.3 


B° - K\ 


1.6 x 10" 


-3 


1.3 x 10" 


-3 


9.3 x 10" 


-4 


0.29 




-32.6 


-38.1 


19.8 








0.41 




0.40 




0.88 




-0.3 


-0.3 


-0.3 


5° -► 






0.37 




0.39 




0.82 




-0.3 


-0.3 


-0.3 


B° -► vrV 


4.4 x 10" 


-5 


4.4 x 10" 


-5 


5.3 x 10" 


-3 


5.3 x 10- 


-3 











5° -► vr !? 


3.3 x 10" 


-5 


3.3 x 10" 


-5 


4.0 x 10" 


-3 


4.0 x 10- 


-3 











B^^rjr) 


5.3 x 10" 


-4 


0.11 




0.11 




0.42 




-1.4 


-1.5 


-0.8 


5° 777^ 


5.9 x 10" 


-5 


0.092 




0.092 




0.13 




0.9 


0.7 


0.5 


5° 7/7/' 


2.2 x 10" 


-5 


0.39 




0.39 




0.93 




0.06 


-0.1 


-0.2 
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TABLE V. Branching ratios in unit of 10~ 5 ,and CP asymmetry in unit of %. 'QCD' and 'EW 
present the QCD penguin and EW penguin effects respectively, and 'DIRAC presents the results 
with the equation of motion. 



AT e ff — 






Di ancnm^ 


^ rvaiio 






CP Asymmetry 


uecay ivioae 


TD 

1 it 






th VV 






QCD 


EW 


DIRAC 


R — , _o_ — 


U.Zo 




U.Zo 


U.Z4 




n ok 
U.ZO 


0.05 


1.9 


4.7 


B~ -»■ 7r~r/ 


U.U40 




U.UoZ 


U.UOZ 




1.0U 


25.2 


0.2 


0.2 




U.14 




U.Z7 


U.Z ( 




l.zO 


29.5 


0.5 


0.9 








n n/i a 
U.U4o 


n n/in 
U.U4y 




n i n 
U.1U 


12.5 


12.3 


11.9 


R — , _0 ZV — 


n m q 
u.uiy 




n i s 
U.lo 


U.Zo 




n i^s 
U.Oo 


-21.1 


-14.9 


-9.5 


R— , JS—J 


U.U4U 




U.4Z 


U.oo 




i kp. 

1.00 


-9.0 


-13.6 


-6.0 


B u ^ K V 


n noo 
U.Uzz 




U.U t\) 


U.U04 




n no/i 
U.UZ4 


10.0 


-8.8 


20.9 


R— , J^O,— — 






n ^4 
U.04 


U.oo 




1.1Z 


-0.3 


-0.3 


-0.3 




U.oO 




l.UU 


i nn 
l.UU 




1 1 o 

l.lz 


13.6 


13.6 


18.7 


r0 _^ ^0^0 
■"d — > 7T 7T 


U.Ul ( 




U.Uo4 


n nQO 
U.Uoz 




n riQn 
U.Uoy 


39.1 


42.0 


46.4 


r0 , _0 ' 


Z.o X 1U 


-3 


n non 
U.UZU 


U.U10 




(1 (11 7 
U.Ul ( 


61.3 


78.0 


78.0 


dO . -.0^, 
B d ^7T 7] 


4.0 X 1U 


-4 


n HQ A 

U.Uo4 


n hq 1 
U.Uol 




n a i 
U.4 ( 


26.1 


29.2 


14.7 


R _^ ' ' 


1.0 X 1U 


-3 


i o v i n~3 
1.5 X 1U 


i o w i n _ 
1.5 X 1U 


-3 


n no(! 
U.UzO 


-6.3 


-7.3 


-6.5 


B d -> vv 


n m q 
U.Ulo 




U.Ul r 


n m q 
U.U18 




n n^n 
U.U (U 


-1.3 


-0.6 


1.1 


r0 , 


U.U < y 




n nan 
U.UoU 


n nan 
U.UoU 




n oi 
U.Zl 


-8.1 


-8.1 


-24.1 


B d ^ K K 






U.UOo 


n nfi/i 

U.U04 




n nan 
u.uyy 


12.3 


12.1 


11.9 


r0 , _+ j^— 

— > 7T A 


U.UOO 




U.ol 


n on 
u.zy 




n kq 
U.Oo 


-28.4 


-30.4 


-19.5 


— > JT A 


i.y x iu 


-3 


U.Z / 


n i c 
U.lo 




n 

U.OZ 


—5 9 


—8 5 


—4.4 


eg - kS 


0.024 




0.44 


0.42 




1.77 


3.1 


0.5 


-0.3 


B° d -► K°r/ 


3.4 x 10" 


-3 


0.027 


6.8 x 10" 


-3 


1.5 x 10~ 3 


-7.8 


-63.9 


-79.6 




0.91 




1.10 


1.10 




1.19 


13.5 


13.6 


18.7 




0.029 




0.055 


0.052 




0.049 


36.6 


39.2 


36.4 


4° - X°r?' 


0.012 




0.021 


0.020 




1.27 


-48.6 


-50.1 


1.7 


5° - K\ 


0.031 




0.037 


0.035 




0.24 


14.2 


9.2 


-16.4 








0.54 


0.56 




1.17 


-0.3 


-0.3 


-0.3 


5° -► 






0.48 


0.45 




0.96 


-0.3 


-0.3 


-0.3 


B° -► TT !?' 


8.7 x 10" 


-4 


8.7 x 10~ 4 


9.3 x 10- 


-4 


9.3 x 10~ 4 











5° -► TT !? 


6.5 x 10" 


-4 


6.5 x 10~ 4 


7.0 x 10- 


-3 


7.0 x 10~ 3 











B^^rjr) 


0.011 




0.10 


0.10 




0.43 


5.1 


5.0 


1.4 


5° 777^ 


1.2 x 10" 


-3 


0.12 


0.12 




0.17 


-4.4 


-4.5 


-3.5 


5° 7/7/' 


4.4 x 10- 


-4 


0.50 


0.50 




1.2 


-1.9 


-2.0 


-1.2 
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FIGURES 



a 




FIG. 1. Diagrams for the matrix elements (P | (7i7 M 75(72 I 0) (fig la) and (P | (717592 I 0) (fig lb). 
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FIG. 2. Leading twist diagrams in QCD for the matrix elements {P\qi'~f^b\B) . 



1 -x)p) P 



1 



-> 3— ► 



4 



q 



P 



A 

B 

\1 



1 ~.x)p 



: 1 



k 



( 1 — y)q 

a 



xp ^ ( 1 — y)q 



FIG. 3. Leading twist diagrams in QCD for the matrix elements (P\qib\B). 
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